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Abstract. Random skew plane partitions of large size distributed according to an ap- 
propriately scaled Schur process develop limit shapes. In the present work we consider 
the limit of large random skew plane partitions where the inner boundary approaches a 
piecewise linear curve with non-lattice slopes, describing the limit shape and the local 
fluctuations in various regions. This analysis is fairly similar to that in OR2 , but we do 
find some new behavior. For instance, the boundary of the limit shape is now a single 
smooth (not algebraic) curve, whereas the boundary in |OR2| is singular. We also observe 
the bead process introduced in [B] appearing in the asymptotics at the top of the limit 
shape. 
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1. Introduction 

A skew plane partition with inner boundary being given by a partition A is any array 
of positive integers tt = {iTij} defined for pairs such that i > 1 and j > Aj, which is 

The third and fourth authors were partially supported by the NSF grant DMS-0354321. 
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non- increasing in both i and j. A skew plane partition is called bounded by N and M if 
TTij = when i > N or j > M. An example of a skew plane partition with inner shape 
A = (3, 2, 2) is given on Figure [2] Looking at this picture we can also see that a bounded 
skew plane partition can be regarded as monotonic piles of cubes in a certain "semi-infinite 
room" with no roof (see also Figure [I]). 

A natural question to ask is: what is the shape of a typical pile containing a large 
number of cubes? One must first make this precise by fixing a probability distribution. 
One natural distribution to choose would be the uniform distribution on all skew plane 
partitions with a fixed number of cubes. This is actually quite difficult to deal with, and we 
instead consider the following probability measure (grand- canonical ensemble) on bounded 
skew plane partitions: 

(1) P(ir) = y^, 

where \ir\ = Ylij n i,j ^ s the total number of boxes in the piles corresponding to the plane 
partition n, < q < 1 and 

TT 

is the normalization factor known as the partition function in statistical mechanics. The 
larger q is, the more likely it is to have many cubes. We will study certain limits of this 
system as q approaches 1, the bounds M and N approach infinity, and the partition A 
defining the back wall of the room grows in a specified way. 

There are many ways to take such a limit. For instance, one can specify that, after an 
appropriate rescaling, A approaches a fixed curve. The limiting behavior of the system will 
certainly depend on this curve. Some such cases have previously been studied: 

• The case when A is empty was studied in the unbounded case in [CK] and in the 
bounded case in [ORlJ. 

• The case when A approaches a piecewise linear curve with lattice slopes was studied 
in [QR2] , 

In this article we consider a limit where A approaches a piecewise linear function with a 
(non-lattice) slope staying strictly between —1 and 1. The large size limit is studied via 
correlation functions using the formalism and notations developed in [OR1, OR2J. 

Before we start, let us mention some other related work. Using results of this article, 
the limit where A approaches a piecewise linear function with arbitrary slopes has been 
studied in [M] . Local fluctuations and the free energy for the partition function on a torus 
were first computed in [NHB . 

1.1. Random skew plane partitions and Correlation functions. Our piles of unit 
cubes are in M 3 with coordinates (x,y,z) centered at points (Z + |, Z + ^, ^ + \) and 
bounded in (x, y) directions as it was described above. That is, the (x, y) coordinates of 
centers of cubes satisfy conditions: 



\\<x<N-\, 

[K+i/2 <y <m. 
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Figure 1. An empty skew plane partition. Here the inner shape is given 
by the partition (6, 3, 3, 3, 1). 



The partition A is given and describes the configuration of the back wall. Denote this 
region of M 3 by D\,m,N- 
The mapping 

/ x + y, 
(x,y,z)^ [t = y-x,h = z — ) 

projects M 3 to M 2 . This projection maps piles of cubes to tilings of the region 

E MiN = {(t,h) | — M < t < N;h> max{-M - t/2, N + t/2}} 

which differ from the one corresponding to the empty pile in finitely many places. This is 
evident from Fig. [2} 

A convenient characterization of a system of random skew plane partitions is the col- 
lection of local correlation functions P(ti,hi),...(t k ,h k ) giving the probability to see horizontal 
rhombi at (ti, hi), . . . , (t k , h k ): 

k 

P(h,hl),-,(tk,h h ) =^2 P (' K )Yl a ia,ja^) 

tt a=l 

where (Jt,h{^) = 1 when (t,h) is the position of the center of a horizontal rhombus and 

o't.hM = ^ not - 

An explicit formula was found in |QRH lOR2] for such correlation functions. To describe 

this formula, we first encode the boundary of the partition A as a function h = —b(t), where 
(2) b(t) = - Y / \t- Vi \--Y l \t-Ui\. 

i=l i=l 
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FIGURE 2. A skew plane partition represented as a stepped surface. Notice 
that the same picture can be viewed as a tiling of a portion of the plane by 
three types of rhombi. 

The points m, . . . , n n _i (resp. v\, . . . , v n ) represent the inner (resp. the outer) corners of 
the back wall. See Figure [T] It is convenient to set no = — M and u n = N. The lines 
t = uq and t = u n are limiting on the left and on the right the region tiled by rhombi. 

Then we introduce the partition of the horizontal interval Z + \ n (no, u n ) into D- U D+ 
where 

D + = {m G Z + - n (no, u n ) \ where the boundary A in (t, h) has slope — - at m}, 

D~ = {m G Z + 1/2 n (no, n n ) | where the boundary A in (t, h) has slope +- at m}. 
Define the functions &±(z,t) by: 



The following holds: 

Theorem 1.1. |OR2] Correlation functions are determinants of the correlation kernel: 



$+(z,t)= H(i-zq m ), <mm)= n (!-^v 



) 




(3) 



P(t u hi),...,(t h ,h k ) = det[K((t i ,/i i ),(t i ,/t i ))] 



l<i,j<k 



The correlation kernel K is given by 




dzdw 



z — w 



zw 
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where C\ is a simple positively oriented contour around such that its interior contains 
none of the zeroes of&+(z, ti) . Similarly C2 is a simple positively oriented contour around 
such that its exterior does not contain zeroes of Q-{w,t2)- Moreover, if t\ < t2, then C\ 
is contained in the interior of C2, and otherwise, C2 is contained in the interior of C\. 

1.2. Main results. We are interested in the behavior of random skew plane partitions 
sampled according to the probability distribution ([I]) for large domains, meaning that a 
limit of the following form is taken: 

• Fix a sequence (r/%) of positive real numbers such that lim^oo = and for all 
k, let q k = exp(-r/ c ). 

• Choose sequences (M&), (iV&) and (Afc) such that as k goes to infinity, the 3- 
dimensional domain rkD\ kj M k ,N k converges to some region D of K 3 . 

• Find characteristic scales of fluctuations for large k and describe random plane 
partitions at appropriate scales in the limit k — > 00. 

The increasing sequence of domains where r& = | and the corners u[ and describing 
the partition are given by: 

uf ] = kUi, vf ] = Wi, 

where n and Uq <V\ < ■ ■ ■ <V n <U n are fixed integers was studied in [0R2J. As k — > 00, 
the function (x,y) h-> hk(x,y) = irkx,ky/k converges in probability to a fixed limit shape. 
This function is infinitely differentiable everywhere except along a curve where the surface 
it defines merges with the walls. Along the curve, its second derivative is discontinuous 
(the first derivative has a square root singularity). The projection of this curve to the 
(t, /i)-plane is called the arctic circle or the boundary of the frozen region. Indeed outside 
this curve the tiling by rhombi is "frozen", which is to say that with probability tending 
to 1, the tiles outside of the curve are of only one type. 

The results in [QR2| are obtained by using the steepest descent method to evaluate the 
integrals in Equation Q. In the limit one obtains exact formulas for the limit shape, as 
well as descriptions of the behavior of local correlation functions. 

In this paper, we focus on a different type of increasing sequences of domains. We keep 
rfc equal to \, but we assume that the partitions A& describing the back wall consist of n 
segments growing in length, and that each of the segments is not straight anymore, but is 
a "staircase" with fixed slope. To be precise, the j th segment of Afc consists of kLj copies 
of the corner with sides (aj,bj) for some Lj,aj,bj G Z>o fixed, see Figure [3j 

In the framework described above, we compute the limit shape of the random pile of 
boxes as k — > 00 and correlation functions. We show that: 

• The system rescaled by a factor r& converges to a deterministic limit shape. This 
is a smooth transcendental function in the bulk. The frozen boundary consists of 
a single smooth curve. 

• As in [OR2], correlation functions in the bulk are given by the incomplete Beta 
kernel. 

• For most values of t, if one fixes t and allows h to approach infinity, the slope of 
the limit shape turns vertical, and the determinantal process with the incomplete 
Beta kernel describing the statistics of horizontal tiles degenerates into the bead 
process introduced in [B]. 
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Figure 3. The setup 



1.3. Relation to Dimer models. The limit shapes for rhombi tilings of bounded do- 
mains with piecewise linear boundary parallel to axes h = const, h db t/2 = const were 
studied in [KOj where it was shown that uniformly distributed tilings develop a limit 
shape given by an algebraic curve. This generalizes the arctic circle theorem for large 
hexagons |CLP1 [J] . Local correlation functions for large hexagonal domains with the uni- 
form distribution of rhombi tilings were studied in (UP, and in [Ke2] . 

The rhombi used in these tilings can be thought of as dimers, since they are the union 
of two adjacent faces of the triangular lattice (monomers). These random tilings are thus 
examples of dimer models. Dimer models have a long history in statistical mechanics, see 
for example |Ka| ITF[ IMC W] . Corresponding tiling models with height functions were first 
studied in the physics literature, where many ideas about the limit shape and the structure 
of correlation functions in the thermodynamical limit can be traced, see for example [NHBJ. 

1.4. Acknowledgments. We are grateful to A. Borodin, A. Okounkov, and P. van Mo- 
erbeke for interesting discussions. 



2. Setup: a piecewise periodic back wall 
We are studying a system of random skew plane partitions, as described in Section [TJ 



The precise system is described in Section 1.2. Let us start by introducing some more 
notation. 

Let t = rt and x = r h be the rescaled coordinates. The r coordinates where the slope 
of the back wall changes are denoted by Ao, A%, . . . , A„, where Ao and A n are the bounds 
of the shape. The t coordinates that correspond to these corners will be A±/r, . . . , A n /r. 
Between A^-i/r and A^/r we have a succession of patterns made of ascending steps 
and bk descending steps, repeated L^/r times. The slope of the piece of wall between Ak-i 
and Ak is = 2 (a^+b k ) • ^ w ^ ^ e convenient to define ao = — \ and a n+ i = \. Let j be 
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the integer in {1, . . . , n} such that r G (Aj-i, -Aj]. Define the function B(t) by 

(5) - J3(t) = "2^0 + ^a/c(^4fc - A fc _i) + oy(r - = - ^^ fc (a fc+ i - a*.) + ayr, 

fe=l k=0 

if < r < Ay B{t) is the limiting rescaled boundary function b(t) of Equation ^ 

expressed in the rescaled variables r and x- We can decide to define B{t) from the 
rightmost corner A n instead of the leftmost one Aq, yielding an alternative expression: 

j n n 

(6) - B(t) = -A n - ^2 a k{Ak ~ A k -x) - (Xj(Aj - t) = ajr + ^(a fc+1 - a k )A k . 

k=j+l k=j 

We are interested in the correlation functions in the scaling limit when lim r _>+o ft\ = 
lim r _5>_|_o rt2 and t\ — t 2 = A(t) is a constant. 

2.1. The function S T;X (z). We want to use the steepest descent method in order to 
find the asymptotics of local correlation functions as the size of the system increases and 
r — > 0. The first step in doing this is to rewrite the integral defining the correlation kernel 
in Theorem 11.11 as 

1 f f <l fel W~<| fe2 W 1 dz dw 



K(( tl , hl ), (h, h 2 )) - {2m? Jr ^ 2i7tw , 

show that s[ r ^{z) converges to some function S TtX (z) as r — >• 0, rt — > r, rh — > x an d find 
the critical points of z h- > S T<x (z). 
The function S^(z) is given by: 

(7) S^(z) = rlo g (z- h - b ^~^ t) 



= -r(h + b(t)) log z + Yl rlog(l-z-V w )- ^ r log(l - ze~ rm ). 

m<t,m£D~ , m > t , m £ 

D+, 

Here and in the sequel, log denotes the branch of the logarithm with an imaginary part in 
(— 7r, 7r), with a cut along R_. 

In order to find the limit S TX (z) of this quantity, we first state a lemma: 

Lemma 2.1. If [C, D] corresponds to a piece of wall with parameters a, b then 
(8) lim V r\og{l-z- 1 e rm ') = -/ log(l - z~ x e v )dv 

C . D JC 
r r 



and 



(9) lim V rlog(l-ze- rm ') = -/ log(l - ze~ v )dv. 

r->0 t—' a + b J (j 



— <m'< — ,m'eD + 

r r 

m'eZ+7; 
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Proof. The proof is elementary. Here is how it goes for (|8|): 



lim rlogfl — z 1 e rm ) = lim rlogfl — z 1 e 



1 „m\ 



r— >0 z — ' r— >0 

C , D , C<m<D,m€D- 
— <m'< — ,m'eD" , i 

r r i m&r(Z+±) 

o-l 

r->0 



lim ^Vlog(l — z 1 e r 



C<m<D, k=0 

mer(a+fe)(Z+|) 
1 a— 1 

= -Vlim V r(a + b) log(l - ^- 1 e m+fcr ) 

fc = C<m<D, 

mer(a+b)(Z+|) 
defn. of/ 1 °" 1 rD .. ,D 



a + 6 

fc=0 



V / log(l - 2~V)cfo = — —7- / log(l - z~ x e v )dv. 
t^Jc a + bj c 



The convergence is uniform in z on compact sets of C \ [e , e ]. □ 



Using Lemma 2.1 for each interval ^4fc + i], [Aj_i,t] 3 [t, and the fact that 

lim rb(t) = B(t), 

rt— >t 

we get 

S TtX (z) :=HmSgW = -(* + 5(r)) log z 
(10) +V — ^— ] g(l-z- x e v )dv + -^- r f T logfl-^VWw 



/ 'logfl -z e -")d w - V fc - / log(l -ze- v )dv. 



0"i ~\~ bn It ^- — ' CLh + bh I A, , 

fc=J-+l K K • /A fc-1 

3. Critical points of 5* TjX (z) 

The localization of the critical points of the function S T x is crucial in order to apply 
the steepest descent method. The following proposition gives a complete understanding of 
their nature. 

Proposition 3.1. Fix Aq < r < A n . Then, 

• there is a unique real number Xo( T ) such that 

(i) For any \ > Xo{ T )> ^r,x( z ) ^ as exactly two non-real critical points, which are 
complex conjugates. 

(ii) For any \ < Xo( r ); S TtX (z) has only real critical points. 

• If t 6 (Aj-\,Aj) ,as x ~~ +°° the two non-real critical points z c ,z c have the 
following asymptotic: 

(11) z c ,z c = e T exp [-e- x p{T)e ±i ^+ a j)^ (i + o(e~ x )), 
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where 



p(r) = n 



k=0 



2 sinh 



T-Al 



Ofc+i-«fc 



and z c is chosen to be the critical point with positive imaginary part. 
If t = Aj + 5, x — > +°° an d 5 — > such that 



x ij) \M 1+a j- a J+i 



is fixed, with 



k=0 



Aj Ah . 
2~ 



Otk+l— a k 



then the pair of complex conjugate solutions behave as z = e T s ^ where s is a 
solution to the equation 



(12) 



The proof is rather technical, and requires a number of intermediate statements. It will 
be completed at the end of this section. 



3.1. An expression for 4zS TtX (z). Using the identities 



a k 1 b k 1 , 1 - ze v , . z - e v , 

ak + bk 2 ak + bk 2 y l-ze u/ x z-e u/ 



and the formula ^ for B(t) , one gets the following expression for the derivative of S 
with respect to z: 
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z-^S(z) = - ( X + B(r)) + y / 1 " <fa + —V / / " dv 

a j + °j i--ze v k ^ +i a k + 6 fc J Akl 1 - ze v 

Ofc ( z ~ e \ a i ( z — e 



f^ak + bk \z- e Ak -^ J Oj + bj 

bi , / 1 — ze _y4 ^ \ 6fc , / 1 — ze 

+ ~ T7T log 1 — + E ~ T7T lo § 

n • 4- hj \ — vp ~ I ^— ' ni 4- n. 



z — e A ^ x 



a i + b, I 1 - ze 7 / <ik + bk \1 — ze Ak - 

fc=jf+i 

i-i 

a fc I z ~ e ' lfc \ a 7 j z — e' 



- (x + B[r)) - £ ^ log ( " ^7T log 

^ afc + h \z- e^-i / aj + bj 

+ Ar^f— V E 

aj + fej \z-e T J + \z-e A k-i 

— xr(^-^)- E -zrto-^i) 

aj + bj , ~ afc + Ofc 



z — e^J- 1 



b(t)+5 

/ \ , t i— m ( z-e Ak \ ,1 , , ( z-e 1 

( 13 ) = " (X " j ) " E(s + bg ( z-e^ J " ( 2 + bg 

1 / z — e Aj \ n 1 ( z — e Ak \ 

+ ( 2 -«,)iog( 7 ^J+ E( 2 -«^^^J- 
> / fc=j+i v y 

Note that each of the terms log ^ Z ~A t *i 1 ) is a holomorphic function of z on C \ 

[e" 4 *- 1 , e Ak ]. Its imaginary part is in (— tt, tt) and represents the angle between the vectors 
z — e" 4 * and z — e" 4 * 1-1 . This angle is when z£R \ [e" 4 * -1 , e Ak ], and approaches tt (resp. 
—tt) when z approaches (e" 4 * -1 , e Ak ) from above (resp. from below). See Figure |4j 

3.2. Asymptotic of critical points when x goes to oo. Here we will prove the last 



z — e A i~ x 



two statements of Proposition 3.1 concerning the asymptotics of the non real critical points 
when x is large. 

Lemma 3.2. Fix Aq < t < A n . Then, for sufficiently large x> S TtX (z) has exactly two 
critical points, which are complex conjugates and have non-zero imaginary parts. Their 



asymptotic as x ^ +°° * s described by the last part of proposition 3. 1 



Proof. First case: Aj_i < r < Aj. 
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Critical points are solutions to the equation 



(14) 



z~S(z) = 0. 
dz 



The real part of this equation can be written as 

3-1 



= If ( z-^S(z] 



-(x- - Xta + a fc) lo £ 



+ (- - ay) log 



fc=l 

z — e" 4 ? 



z — e 



z — e* 



'2 + a i) lo £ 



z — e 



2 — e^- 1 



z — e' 



n 1 



fc=j+l 



z — e 



z — e 



(X - \) ~ logjl - 2e r | + ^Ofe+i - a fc )log|l 



ze 



-A fc | 



k=0 



From here we can see that a solution z as x ~ > +°° should approach to e T , or to e Ai , 
for some I £ {1, . . . , n} with < a; (inner corner). 

The imaginary part of equation ( 14 ) for the critical points is: 

3-1 



, z — e 



1 1 

(15) = - + a fc )angle(z - e^- 1 , z - e Ak ) - (- + a i )angle(z - e^-\ 

k=0 

1 ™ 1 

+ (- - a i )angle(z - e T ,z - e A i) + ^ (- - a fc )angle(z - e At -\z- e Afc ). 



fc=j+i 



A_i A, 



Figure 4. Geometric representation of the imaginary part of log ^ 

as the angle between z — e Al and z — e 

Assume that z is a critical point with positive imaginary part, and z — > e Al for some 
I (since the equation has real coefficients, the non real solutions will come as pairs of 
conjugate complex numbers). Since z approaches to the real axis, the angle between 
z — e B and z — e c approaches to zero, unless B or C is equal to A\. Therefore as x ~ ► +°° 



, two terms in Equation (15) will dominate. These are 

J (| + «/)angle(z - e^- 1 , z - e Al ) + (~ + a i+ i)angle(z - e Ai , z - e^ 1 ) if r < A;, 
I (| — a/)angle(z - e^- 1 , z - e Ai ) + (| - a; + i)angle(z - e Ai , z - e Ai+1 ) if r > A;. 
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Thus, these expressions should vanish as x ~ * +°o- 
But this is impossible, because the sum 



angle(z - e Al ^ , z - e A[ ) + angle(z - e Al , z - e Al+1 ) = angle(z - e Al ^ , z - e Al+1 ) 

is positive and close to tt (see Figure [4]). Therefore this option for a critical point is 
impossible. 

The only possibility left is that z — > e T as x ~ * +00. Let z = e T ~ 6 , where e — > as 
X — > 00 and arge E (— n, vr). 
As e — > 0, we have for k 7^ j: 



log 



z — e 



fop 



z — e 



z — e 
z — e^- 1 



log 



+ 0(e) 



A k — A k _i 



+ log 2 sinh 



T-Ak 



t-A 



3-1 



+ log(— e) — log 2 sinh 



— log 2 sinh 

T - Aj-! 



+ 0(e) 



log 



z - e A i 



z — e' 



A - — T 

— log(e) + log 2 sinh 



r-A j 



+ 0(e). 



Substitute these expressions in the equation for critical points: 
j-i 



= -(X - S " X>* + i) o^" 1 + ^g 2 sinh 

fc=i ^ 



2 



log 2 sinh 



T ~ Afe-l 

2 



+ (2 ~ a i) 



1 /t — A-—i 

— — h log(— e) — log 2 sinh 



r- A 



— T 

— h log 2 sinh 



r - An 



A k - A 



k-l 



+ log 2 sinh 



- log(e) 
r-A k 



log 2 sinh 



t-A 



k-l 



(16) 



1 ™ T — A 

-X ~ log(e) =F ("j + ^)wr + 5Z( a fe+i - a/c)log2sinh — h 0(e). 



fc=0 



+ 0(e) 



Here we used the identity 



log(-e) 



log e — in if arg e > 0, 
loge + m if arge < 0, 



which holds for our choice of branch for log. The cancellation of the linear combination of 
the Ai comes from the two possible expressions for —B{t) given by Equations §5§ and Q. 
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Solving this equation for e, we arrive to the asymptotical formula: 



n 

fc=0 



2 sinh 



t-Ai 



[l + 0(e-x)) 



P(r) 



which implies that the asymptotic of the pair of complex conjugate solutions is given by 



the formulae in the last part of proposition 3.1 



Second case: r becomes close to Aj for 1 < j < N — 1 as x g° es to +oo. 
When r = Aj + 5, 5 — > and x ~> +oo, the same arguments as above show that complex 
critical points still accumulate near z = e T but instead of the Equation (16) for e, where 



we will have 



(17) X-X 



(f) 



-(ay+ *)log(-e) + (~-ay)log 



(--a i+ i)Iog(e-<y) + 0(e). 



From here we obtain equation (12) for s = e/\5\ 



This completes the proof of Lemma |3.2 



□ 



Whereas the previous lemma describes what happens when x goes to +oo, the situation 
when x goes to — oo is given in the following: 



Lemma 3.3. For each Aq 
S(z) is real. 



< t < A n , and sufficiently negative x, every critical point of 



3.2 



To have ^(z£S(z)) 



0, z needs 



Proof. This is proven in the same way as Lemma 

to be very close to e Al , for some I € {0,1, ... ,n}, such that ati + \ > a\. Exactly as in the 
case x ~^ °°> one can show that if I ^ or n, then *^s(z-^S{z)) cannot be 0. If z is near 



e A ° or e An , then z is actually real and respectively z < e A ° or z > e Ar ' 
when x is sufficiently negative, z4-S(z) = has no complex solutions. 



This shows that 
□ 



3.3. Proof of Proposition 3.1[ Note that the last two parts of the proposition were 
proven in the previous section. To prove the first part we first show that for fixed r G 
(Ao,A n ) there is a unique value of x = Xo{ T ) such that the number of critical points of 
S(z) with non-zero imaginary part changes as x passes through xo- That is: 

Proposition 3.4. For each t 6 (Ao,A n ) there exists xo such that the number of critical 
points of S(z) with ^s(z) ^ is the same for all x > Xo> an d it decreases by 2 when x 
passes the point xo into the region x < Xo- 

We prove Proposition [3^4] at the end of this section, after setting up some technical tools. 



Define f{z) 



' dz 



S{z). It is clear that z = is not a critical point of S(z) and therefore 



critical points are zeroes of f(z) in C \ {0}. It is also clear that z is a double critical point 
if and only if 

/(*) = f'{z) = 0. 

The following two lemmas give some information about the location of the double critical 
points of S(z). 

Lemma 3.5. There are no solutions to f(z) = in [e A °,e An \. 



14 C. BOUTILLIER, S. MKRTCHYAN, N. RESHETIKHIN, AND P. TINGLEY 

Proof. Consider z = x + ie, where x £ {e Ai - x , e Ai ) and |e| << 1. Looking at the terms in 
f{z) in formula (13), it is easy to see that for such z we have 

9 (/(*)- clog (^^))=0(e), 



where c / is the coefficient of log y z ~l i l 1 j m fi z )- Since log is the branch of 

the logarithm with an imaginary part in (— tt,tt) with a cut along K_, we have that 

9 ^clog ^ Z < aI-x ^ = ^ zC7r + 0(e), which in turn implies Q(f(z)) = ±cir + O(e). This 

shows that there are no solutions to f(z) = in (e Ao , e An ). 

The only remaining points are e A ^ but the function f(z) is singular at these points, 
and therefore they cannot be solutions either. □ 

Lemma 3.6. f'(z) = has exactly one real solution outside the interval [e A °, e An ] (possibly 
at oo ). 

Proof. Recall that, by definition, the equation f'(z) = means 



(18) f{z) = { *>__ . d 



d f(z) = z ^ oo 

lim^oc z 2 j z f(z) = z = oo. 



Thus solving f'(z) = in E U {oo}\[e A °, e An ] is equivalent to solving G(z) = in this 
region, where 

(19) G(z):=(z-e A ^)(z-e T )f'(z). 

When r G (Aj-x,Aj) we have 

(20) 



2 +a V ie 6 j ( Z - e r)( z _ e A J _ l) + ^ 2 e \ z _ e A 3){z _ eT) 

+ 4^l 2 °V (e } (z-e^)( 2 - e ^)- 
Notice that each term has the form 

(z - ui)(z - u 2 ) 

Q 

(z - v\)(z - v 2 ) 

where 

C < 0, til, U 2 > V\,V2, OR C > 0, Ui, U2 < V\,V2, 

and 

> til, U2,V\,V2, OR z < Ui,U 2 , Vi,V 2 - 

Under these conditions 

d ( , /(z-in)(2-it 2 )\\ / 1 1 1 1 \ 

— clog ) ^ 2 -± = c + < 0, 

dz \ \{z - Vl){Z - v 2 ) J J \z-ui z-u 2 z-vi z-v 2 J 
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for z G R\[e °, e ™], so G(z) is decreasing in z for all z E M\[e j4 °, e" 4 ™]. Furthermore, it is 
clear from formula ( |20| ) that G(z) is regular at oo. Thus G(z) is decreasing as one moves 
along (e An ,oo], and then "wraps around infinity" and moves from — oo to e A °. 
The function G{z) can also be written as 

/ i i i 

„rW 2+°! a fc ~ Qfc+1 , 2 



G(z) = (z-e A J)(z-e T ) 



-- V 



+ 



Or 



1 



fc=l 



Using this formula and the fact that ^ + a\ > and | — ctn > 0, one can see that 

lim G(z) < 0, lim G(z) > 0. 

this completes the proof of the Lemma. 



□ 



infinity. L emm a 
Lemma 



3.6 



Proof of Proposition 3.4 As x varies, the number of non-real critical points can change 
only at the poi nt xo when there is a double solution to f(z) = on the real line, or at 
shows that this never happens with z G [e A °, e An ]. 
shows that f'{z) = has exactly one solution z$ E M\[e j4 °, e An ] (or possibly 
at oo). Since f(z) = g(z) + x, where g(z) does not depend on x-, it implies that there 
is exactly one point xo = ~g( z o) where f(z) = f'(z) = and that this happens when 
z = zq. □ 



Proof of Proposition ^. l\ By Lemma 3.2 for sufficiently large x there are exactly two 
critical points of S XyT (z) with 7^ 0. By Lemma 3.3, for sufficiently negative x there 



are no such critical points. Taking into account Proposition 3.4, this implies there is a 
unique value xo such that all critical points of S(z) are real for x < Xo and there is a pair 
of complex conjugate critical points z c ,z c with 9?(z c ) > for % > Xo- This completes the 
proof of Proposition 3.1 □ 



3.4. The asymptotic of critical points near corners. The following proposition de- 
scribes what happens in situations interpolating the two regimes described in Proposition 

EH 



Proposition 3.7. Complex conjugate solutions to (12) have the following asymptotic when 
p ->■ 0: 



s = p 1+a j a i +1 exp 



,^ + °'' )lHOM). 

1 + Oij — Otj+l I 



If p —> 00 the asymptotic is different for different signs of 5: 
(i) When 5 -> +0 



s = p- 1 exp r±*7r(i + a j+ i)j (l + 



(ii) When 5^-0 



s = p' 1 exp (±™(\ + otjYj U + {~)) 
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Proof. We will do the case when p — > oo, 5 > 0. The other cases are similar. From ( 12 ) we 
see that if p — > oo then either s — > or, if (<x/+i — otj) < 0, s — > 1. If s — ¥ 0, then s can be 
easily calculated to have the asymptotics given in the statement of the proposition. We 
need to show that the case s — > 1 is not possible. 

We have z = e T ~ £ , e — > 0, and we can assume that Im{e) G (0, tt). 

Write s as s = 1 + te %e . Since Im(e) G (0, 7r), we can assume 9 G (0, 7r) as well. 



Let's rewrite (17) as 
X ~ X U) + (1 + uj - a j+1 ) \og(5) = 

= "("; + \) log(-a) + {\~ <*j) log(^-) - (- - atj+x) log(a - 1) + 0(<5) 

= -(a, + 1) log(-l - ie ie ) + {\- a,) Mj^^) - - log(ie ie ) + 0(<5). 

Looking at the real part of this equation we get 

X - X ij) + (1 + (*j ~ Uj+l) log(5) = - ay) log(i) - (i - a j+ i) log(t) + O(l). 
From here we get 

t = e -)■ 0. 

Let's look at the imaginary part. We have 

= -(a, + ~)(-tt + 0(t)) + (i - - O(i)) - (i - a i+1 )0 + 0(5). 
From here 

or 4- i 

fl = 7T J 2 +0(t)+0(g). 
But this is impossible because 

9 £ (0,vr) 

and 

1 

< aj -a j+1 < dj + -. 

□ 

These asymptotics for s as a function of p imply the following asymptotical formulae 
for complex conjugate critical points z = e T_e of S(z). 

• When x —> +°° an d 5 such that e x |<5| 1+aj-aj+1 — > 0, we have 

i ■ a i + h 
±m-. 1 — ±- 



e = e 1+a i- Q j+i e 1+Q j- Q i+i (1 + o(l)). 
When x ~^ +°° an d ^ ~~ > +0, such that p — > oo, 

e = e ± ^ + 5) e -x|<5p-+ 1 - Q J(l + o(l)). 
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Figure 5. The frozen boundary when the parameters are 
A = {-8, -7, 0, 4, 6, 7}, a = {-0.1, 0.4, -0.45, 0.4, -0.25} 



• When x —> +°o and 5 — > —0, such that p — > oo, 

e = e ±iw ^+ 1+ s) e - x |<y|^+ 1 - a i(l + o(l)). 
Notice that the limit x ~ * 00 °f complex conjugate critical points for fixed r given by 



( 11 ) agrees with the last two asymptotics when r — > Aj ± 



4. Correlation functions, the frozen boundary, and the limit shape 

In this section we will study correlation functions in the limit of the infinitely large 
system. In particular, the one-point correlation function gives the macroscopic density of 
horizontal rhombi. 

In our analysis we will follow [ORl| and [OR2]. The basic idea is to use the steepest 
descent method for computing the asymptotical behavior of the double integral defining 
the correlation kernel. 

It has been shown in [ORl] that if (r, x) is such that all critical points of S(z) are real, 
the region in the vicinity of this point is frozen, i.e. in terms of tiling by rhombi it is 
tiled with probability 1 by rhombi of one type (tilted to the left, tilted to the right, or 
horizontal). 

The region where two simple real critical points collapse into one degenerate critical 
point is a curve in the (r, x)-plane which separates the frozen region from the disordered 
region, where the function S(z) has a pair of complex conjugate simple critical points. 
This curve is the analogue of the arctic circle for tilings of large regular hexagons fCLPj . 



As it follows from Proposition 3.1, and in opposition to what occurs in [DRlJ, the 
boundary of the frozen region in our case projects bijectively to the interval (Ao,A n ) on 
the r axis. As we will see below, the region above this curve is disordered, the region below 
is frozen. 
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4.1. Local correlation functions in the bulk of the disordered region. In this 
section we will prove the following theorem. 

Theorem 4.1. In the limit when r goes to 0, the correlation functions of the system near 
a point (x, t) in the bulk are given by determinants of the incomplete beta kernel 



where the integration contour connects the two non-real critical points of S T)X (z), passing 
through the real line in the interval (0, 1) if A(i) > and through (— oo,0) otherwise. 

The proof is completely parallel to the similar statements from |ORl| IOR2] . Here is the 
outline. 

Proof. The correlation functions are given by Formula Q. When r — > 0, the leading 

asymptotics of the integrand is given by e r T^w- ^ n ^ e region where S T>x (z) has 
exactly two complex conjugate critical points we calculate the asymptotic of the integral 
by deforming the contours of integration as follows. 

The poles of $~(^] are real and lie in the interval (e rtl , oo). The poles of are 

also real and lie in the interval (0, e r * 2 ). When r — > 0, they accumulate along the intervals 
(e T ,oo) and (0,e r ) respectively. 

It is clear that any deformation of the contours of integration away from the real line 
doesn't change the integral, as long as the point where the z-contour crosses the positive 
real line does not move to the right, the point where the w-contour intersect the real line 
does not move to the left and the contours do not cross one another. 

Note that whether the z-contour is inside the u>-contour or vice-versa, depends on the 
sign of At = t2 — t\. Now look at the level curves of 3t(S(z)) passing through z c and z c . 
In the case when the z-contour is outside, we can deform the original contours C z , C w , to 
C' z , C' w as in Figure [6j The only poles we cross are the poles at z = w, so we need to take 
into account the contribution from the residues at z = w. We obtain 



where C is an arc connecting points z c and z c which crosses the real line on the positive 
side. In the case when the t/;-contour is inside of the z-contour, C crosses the real line on 
the negative side. 

Now, let's show that in the limit r — > the first integral is 0. Since lini^o ^tS(z) = +oo, 
the gradient of the real part of S(z) looks as in Picture [6j which implies that everywhere 
along the contours C' Z ,C' W except when z = w = critical points, we have %t{S(z)) < 

$l(S(w)). Because the asymptotically leading term of the integrand is e r i~^> 
we conclude that the integral is as r — > exponentially fast. Since z— and w— con- 
tours intersect transversally, there are no problems with the integrability of in the 
neighborhood of the critical points. 



(21) 




(22) 
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Re(S(z))=Re(S(z cr )) 



Figure 6. Deformation of contours when the z contour is outside the w contour. 



We conclude that the correlation kernel is equal to the second integral in (22). 



The second integral in (22), i.e. 



$_ (z, ti)$+Q, t 2 ) zh2 _ hl+B(t2) _ B(tl) dz 
i e $ + (z,ti)$-(z,t 2 ) 2mz' 
is asymptotic, as r — > 0, Ur — > t, ti — t 2 = At, h\ — h 2 = Ah to: 



(23) 



(-e T ) 



iv_(i0-iv_(t 2 ) 



/ (1 - e-z) 



At z -Ah+f 



t dz 



2mz 



where N^(ti) is the number of descending unit pieces of back wall between the leftmost 
corner Aq and ti . The prefactor in front of the integral is of the form 

g(h,hi) 

g(t2,h 2 )' 

which disappears by multi-linearity when we compute correlation functions, as determi- 
nants of the form det(K((ti, hi), (tj, hj))). 

□ 

The kernel T has a simple expression on the diagonal: when At = Ah = 0, we have 



4, T (0,0)= f 



dz 9 
e 2inz 7r' 



where 9 is the argument of the critical point z c in the upper-half plane. This quantity 
is exactly the density of horizontal tiles. It is therefore the vertical gradient of the limit 
shape height function. 
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The correlation kernel depends only on one complex parameter z c . It is determined 
uniquely by coordinates (t,%) or by the slope of the limit shape at this point. 

Corollary 4.1. The distribution of the tiles in the neighborhood of a point (x, t) converges 
to the translation invariant ergodic Gibbs measure on tilings of the plane with rhombi, with 
activities 1, e~ T \z c \ and \z c \ for the three types of tiles in absence of a magnetic field. 

Indeed, from [KellKOS] the correlation functions for the translation invariant probability 
measure with activities a, b, c on rhombus configurations (with no external field) are given 
by determinants of the correlation kernel 



z -Ah+f w At dz dw 
T 2 a + ^(b + - ) 2inz 2inw' 



Taking the proposed choice for the activities, and computing the integral over w by residues 
leads to a simple integral along a sector of the unit circle. After a change of variable in 



the integral by multiplying by \z c \, one gets the same expression as in (21), see jpRlJ. 
The convergence of finite dimensional distributions follows from Theorem 4.1 Since the 
space of tilings is compact for the product topology, the convergence of finite dimensional 
distributions implies the convergence of the whole distribution. 

This result gives a precise description of the local behavior of the system in the neigh- 
borhood of a point inside the liquid region of the limit shape. However, the expressions 
of the critical points z c and ~z^ are not explicit for a finite %■ I n the next subsection, we 
investigate the asymptotic when \ ~ * 00 • 

Similar arguments can be applied to the case when all critical points are real. In this 
case the integral Q tends either to zero or to one depending on the sign of the critical 
points. The steepest descent method results in the following asymptotic: 

where a > 0, p = 1 when (r, x) is below the segment of the boundary of the limit shape 
between two points where it touches the boundary, and p = when (r, x) is below the 
boundary of the limit shape but above the turning points. 



4.2. The correlation functions for large \ and the bead model. In this section, 
we investigate the behavior of the high piles of cubes of the random skew plane partition, 
that are close to the back wall. The asymptotic of local correlation functions as x ~ > 00 
depends on whether r is in a vicinity of Aj or not. 

4.3. Aj-i < t < Aj. In this subsection we consider the case when A^\ < r < Aj. 
Let 

ti = - + T]i, hi = -+e x ii, £ = 1,2, 
r r 

so that Arj = At and A£ = e~ x Ah. 



Theorem 4.2. In the limit where x ~^ °°; such that Ar] and A£ are fixed, the kernel 
K®. (£2,^2)) of Equation (21) divided by e x and a factor which does not affect the 
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determinant M), converges to 



^ (7) ((m,ei),(%,e 2 )) 



V1-V2 



dip 
2^' 



V1-V2 



-«<M5i-6) 



dip 
2^ 



m > m, 



where ctj is the slope of the piece of the back wall between Aj-\ and Aj, 7 = sin-Tray, and 



d = p(r) cos(iraj) = cos(7raj) JJ 



k=0 



2 sinh 



t- Al 



Proof. Suppose t\ > £2- Recall that the two complex conjugate critical points z c and z c 
for this value of r have the asymptotics (11) when x — >• 00. The curve C joining z c and z c 
can be chosen as the positively oriented arc of the circle centered at 0, of radius \z c \. A 
possible parametrization of this arc is: 



z = z 



(ip) = e T exp (e x p(r) (sin(7raj) + iip cos(naj))) (1 + o(l)), 



where ip runs from —1 to +1. 



The expression of the kernel (21) becomes 



K h T (At, Ah) = e' x p(r) cos (7^ 



'1 



- , >\Ai/ , ^-Ah+At/2dip 
(*(¥>)) 



When x goes to infinity, one has 

(1 - e~ T z{ip)) ll ~ t2 = (-e-V(r)) ,)1 "' e (tan(xaj) + i<p cos( TOj )) r ' 1 "' 72 (1 + o(l)), 

The factors 



(e~ x p(T)) m ~ V2 ^e TeX e p ^ sin ^ a ^ 

are of the form . These factors have no effect on the computations of the probabil- 

ities, since they cancel out when computing, in the limit, the determinant det(i^ X)T (tj — 
tj,h{ — hj)). We can thus drop them out, and obtain a new kernel defining the same 
determinantal process. 

The integral converges as % ~ t +°° to 

p{r) cos^o,) C (sin(7ra j ) + iip cos^ay))™ e -ifP(r) oo^a^-^) 



2vr 

Similarly, when t\ < t2, we use the same change of variable z = z(ip). But now, 



/\\..(A/. A//i divided by ^77-77; 
— p(r) cos(7raj) 



t. ' i - s 1. — y)- The same computations as above show that 

p(r) cos(7rcrj) ' / V ' p(T) cos(iraj) ' ^ 

'- converges to 



-1,1] 



(sin(7r<x,-) + iv9cos(vra i ))'' i e"'=''' 1 ' 



»H-»B r -i<pp(r) cos(7ra,)(6-g2) ^ 

2tt' 
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The kernel divided by e~ x and some factors 9 \ Tll 'f 1 \ which cancel in the determinant, 
converges to 

P (t) cos(7rQ,-) f (sin(7ra,-) + i<p cos^a,-))" 1- " 2 e - <CM (™>( T K6-&)*'^ J tl > t 2 

J-l 27T 

-p(r) cos^) / (sin^a,) + icp cos(7ra j )) ,?1 - r ' 2 e -*co<™y)/>M(?i-6)<^ ^ < fe 
7r\[-1,1] 27 i" 

This ends the proof. □ 

As a consequence, using the same argument as in [B], we have the following theorem 

Theorem 4.3. The point process describing the positions of horizontal rhombi in the 
neighborhood of a non frozen point (r, x), with < r < Aj in the limit shape r — > 

converges to the bead process onZxM with parameter 7 = sin7raj and density 



1 1 / — A \ 

— p(r) cos(7ray) = — cos(7ra,) I I 2sinh I — - — I 

7T 7T - LJ - V 2 / 

fc=0 



Note that the density is singular when at the corners: it tends to at inner corners 
(where otk+i > ctk), and goes to 00 at outer corners. This is a remnant of the singularities 
of the limit shape observed in |ORlj . This suggests that the scaling has to be modified to 
observe a non trivial phenomenon in the vicinity of the corners Aj.. 

4.3.1. r is in the vicinity of Aj for < j < n. Now assume that r = Aj + 5, and 5 — > 
as x ~ * 00 i n such a way that p = e^^' 1 \5\ 1+aj ~ aj+1 is fixed. Recall that in this case, the 



critical point of S T , x (z) are given by Equation (12). In this limit the curve C joining two 
complex conjugate critical points can be parameterized as 

z(4>) = e r exp(-^e-^ ( *-* (j)) (s' + is" (/>)), 
where s' and s" are real and imaginary parts of s respectively, and /3 



J l+atj—ctj+i ' 

Theorem 4.4. Assume that x — > 00 and 5 — > as above, and that An = At and A£ = 
e~P jX Ah are fixed. Then the correlation kernel (21), divided by e _ ^ x and some factors 



not affecting the determinant, converges to the same expression as in the previous theorem 

, 171—772 



with 



'd! f\ (y + 1 2 e-^(fr-&>^ if vi > m 



d' = e^ U) s" 



s" 



s 



The proof is completely parallel to the proof of the previous statement, so we will skip 

it. 

It is easy to see that these formulae when 5 — > ±0 and p — > 00 agree with the previous 
result when r — > Aj±0. let us verify it for 5 — > — 0, the calculations being almost identical 
for 5 —> +0. We have in the limit 5 —> 0—: 



s ~ p e 



l p i7r(Qj+i) 
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From the definition of p and a simple algebra we have 



(pe-x)ft-i = |j| ^ e - 
Taking into account the definition of A£ and A£ we obtain: 

ACe^s" ~ |$|«J+i- a » e ft* W cos(vra i )Ae 



Moreover, when r — > A j — 



p(r) ~ ^)|(5| a i+i-"i = |<y| a i+i- a i e * tf) . 
where 5 = t — Aj. Comparing the two kernels in these regimes yields: 

e-ft^ , )((7 7 i ) Ci),(r?2,C 2 ))~e-^W((r ?1 ^ 1 ),(r ?2 ,6)). 
It is also easy to find the asymptotic when t — > 0. In this case: 

/ • ( 77 + a j 
7 ->sm — — 

\2 1 + aj — aj+i 

and for the density we have: 

M 3) t^s" cos(^ a ^ + a i ). 

2 1 + ay — q j+ i 



Again, as a corollary of Theorem 4.4 the statistics of the horizontal tiles, rescaled 
vertically by a factor e _x ^ , converge as x ~ ► 00 to the bead process with parameter 

7 = - ( a ff£lffi) ) and density J (pO))^ cos (^^)) ■ 

Notice that in the vicinity of r = one can get a bead process with any parameter in 
the interval [sin(7ray), sin(7raj+i)] by appropriately tuning p. 
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